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Abstract 

Using the statistical description common to random waves and quantum fields we show how the 
probability of having a nodal line close to a (translationally symmetric) reference curve 7 is related 



to the Casimir energy of an appropriate configuration of conductors. 
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I. INTRODUCTION 



Random waves (RW) have been an object of interest for their statistical properties both 
in wave mechanics and statistical mechanics. In wave mechanics they turned out to be an 
incredibly interesting and rich scenario for studying the statistics of topological properties, 
ike phase singularities |l] . In optics they made a good statistical model for speckle patterns 
3] in laser beams. In quantum mechanics they have been studied P, |^ in connection with 
semiclassical wave functions in chaotic billiards. In statistical mechanics their properties 
have been put in connection with the statistic of defects and vortices nf and an interesting 

n 

duality with a percolation problem has been put forward recently p. 

The purpose of this paper is to exploit in a new direction the description of RW in terms 
of quantum field theory (QFT). 1 I will make a quantitative connection between the ground 
state energy (or Casimir energy) of a scalar field in a given configuration of semi-penetrable 
conductors in d — 1 space dimensions and the probability of having a certain configuration 
of nodal lines in d dimensions. 

This paper is far from being exhaustive or self-contained. I will briefly introduce the 
concept of random wave referring the reader to the existing literature for their statistical 
properties, in particular the properties of their nodal lines. I will then rephrase the concept 
of Casimir energy for a scalar field in a static background (a toy model for QED where 
the static background models the conductors) in a language closer to that of random waves 
statistics. I will then point out the connection between Casimir energy in this background 
and probability of having a certain nodal line configuration. Finally I will draw some con- 
sequences on the nodal line probability and comment on the possible extensions on which 
further work is needed. 



II. RANDOM WAVES 



An isotropic random wave (RW) in d dimensions is the random function defined on a 
subset of M. d (we will not use any particular notation for vectors but there is little room for 



The statistical mechanics description of RW is related to this via the usual QFT-statistical mechanics 
duality (Wick rotation). In this sense this description is already contained in the work of B. Halperin in 
0] who employed it to study the statistic of vortices and defects. 
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confusion) as 

J [2 

4>{x) = J ~j e (kj) cos(kjX + 8j) (1) 
i=i 

where the phases Sj are uniformly distributed in [0, 2tt) and the vectors kj are random 
variables as well. We will assume isotropy of e, i.e. e(k) is an even, analytic function of the 
length of the vector k. 

For any finite J the moments (0(xi)...0(x n )) are not factorizable, but in the limit J — > 00 
Wick theorem holds [l| (among other things one requires the existence and finiteness of at 
least the second moment, i.e. (0 2 (x)) < 00): 

(<P{x 1 )..4{x 2n )) = (^ x M x i)) - (<P{xk)<P{xi)) . (2) 

Contractions 

In the following we will hence always assume the limit J —>■ 00 is taken. 

Wick's theorem is equivalent to saying that the statistical properties of RW can be de- 
scribed by a Gaussian probability functional 

P[(j)] = I ex p J d d xd d x'<P(x')h(x' ,x)<P(x)^J , (3) 

where Z is a normalization constant and h(x',x) = h(\x' — x\) for isotropic RW. From this 
probability functional the reader could already recognize the usual set-up of the statistical 
mechanics of a non-interacting real field 0. The function h is determined by the spectrum 
e(k) (and vice versa). We will now determine their connection. 

To this purpose is convenient to pass to the Fourier components of the field <pk = 
j d d xe lkx 4>(x) and define h(k) through 



d d xd d x'e ikx - ik ' x 'h{x',x) = {27c) d 5 {d \k' - k)h(k). (4) 
In terms of <fik the probability functional is 

^]=^-p(4/(^^i^i 2 )- ^ 

The limiting Gaussian probability functional (jUJ) or (0) can describes the statistical properties 
of (JH) if we choose the spectrum e(k) as 

lim -e 2 (k) = t^-t^t (6) 

J-00 J v ; h(k) (2vr) d V ; 
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which means that in the limit J — > oo the sum over kj must be substituted by the integral 
in d d k whose measure is given by right-hand side of (|B|). This is the promised connection 
between h(x) and e(k). In this way when J — > oo the propagator G tends to 

G(x,0) = (0(x)0(O)) = Jim gle'^Ocoste*) = | -^-^e* (7) 

where we used the fact that e is even in k to substitute e lkx for cos(/ca;). 
There are at least two 'natural' choices for the spectrum h(k): 

• The scalar field spectrum l/h(k) = 8(A — \k\)/(k 2 + m 2 ) where one has to introduce 
the cutoff A to ensure the finiteness of G(x, x) = (4> 2 (x)). 

• The very singular monochromatic spectrum h(k), such that l/(2irh(k)) = S(\k\ — K). 
This last choice gives G(x, 0) = Jq(Kx) which is a statistical model for the solutions 
of the Schrodinger equation — A^ = K 2 ip in chaotic billiards. 

III. CASIMIR ENERGY AND NODAL LINES 

We now turn to the main point of this paper: the connection between nodal lines proper- 
ties and Casimir energy. For simplicity at the moment we assume d = 2, the generalization 
to other d will^be straightforward. 

we introduce the functional 



Following 



X,[<P] = l - I dsrr[ (Si 
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where the integral is defined over the reference line 7 = {x(s)|s e[0, £]} and parameterized 
with the length of the line itself, s. For any given reference curve 7, Xy[</>] is a random 
variable whose generating function 5 7 (A) is defined as 

5 7 (A) = (e~ xx ^) = J V<p P[(j>] e -^ x ^ ds ^ 2( * (s)) . (9) 

It has been shown in 0| that S 7 (X) can be interpreted approximately as the probability of 
having a nodal line in the tube of radius r = (((V0) 2 ) A) -1 / 3 built around the reference 
curve 7 (in d dimensions 1/3 gets substitutes by l/(d + 1)). Notice that the radius r of 
the tube goes to zero when A — ► 00. The approximation allowing us to interpret as 
the probability of having a nodal line relies mainly on a mean-field approximation where 



4> 2 J '(V '4> 2 ) — > 4> 2 / ((V0) 2 ) as discussed in j4J. It is not easy to estabilish the limits of this 
approximation so we will adopt it as a working hypothesis and we will see later a situation 
in which it possibly fails. From now on we will simply say that 5* 7 is 'the probability to have 
a nodal curve 7' without referring to the tube radius r or the approximation within which 
this interpretation has been derived. 

Let us write in Q P[<f>] explicitly inside the probability functional 

S 7 (A) = J Zty-^exp f~ J d 2 xd 2 x' <f)(x')(h(x',x) + 5 {2 \x' - x)V(x))4>(x)j , (10) 

where we have defined 

V{x) = \ fds 5 {2 \x-x(s)). (11) 

J-y 

Let us now specialize the problem in two ways: 

• Choose h(x', x) to mimic a scalar field, with a cutoff A intended in all the momentum 
integrals 

h{x',x) = 5 {2 \x' -x)(-A + m 2 ). (12) 

• Consider a random wave in the strip [0, T] x M. Denote the two cartesian coordinates 
in the plane as Xq, x\ so < Xq < T and X\ G K. Choose the reference line 7 as made 
of n > 1 disconnected lines parallel to the Xq axis and intersecting the X\ axis at the 
points {ai, a n } 

7 = {a u ...,a n } x [0,T]. (13) 
With this assumptions the final expression for the generating function 5* 7 (A) is then 

S y (\) = I / P0exp ( -- f d 2 x 0(x)(-A + m 2 + V(x))<p(x)\ (14) 

where A = d 2 / 'dx^ + d 2 / 'dx\. This expression itself is reminiscent of two intertwined concepts 
in QFT and statistical field theory: the Casimir energy £ in the first and the free energy 
F in the second. The connection with the latter is evident, without any need for formal 
manipulations, F = — log(S 7 (A)). The connection with the Casimir energy becomes evident 
as well if we perform a clockwise (inverse) Wick rotation in the x$ coordinate, xq — > it. Then 
(114)1 becomes 

\ I £>0exp(4 / dtdxi 0(t,Xi)(-<9 2 - m 2 - V(x 1 ))<f>(t,x 1 ) J = e"^ T . (15) 

^ '[0,T]xR / 
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Here d 2 = d 2 /dt 2 — d 2 jdx\ and £ 7 the Casimir energy in the background V and 

n 

V(y) = J2^-y)- (16) 

8=1 

We can now establish the promised connection between the generating functional 5^ (A) 
and the Casimir energy £ 1 of the corresponding background as 

5 7 (A) = e~ T£ -. (17) 

In words: The probability of having a (translationally symmetric) nodal line 7 in a random 
wave ensemble is related to the Casimir energy of a configuration of conductors given by a 
constant-time section of'-f. 

The generalization to d dimension is easily obtained (and is already understood in the 
previous paragraph). Since a nodal hypersurface has codimension 1 so does its constant-^o 
section. Hence the problem maps to the Casimir energy of codimension 1 surfaces. The 
dual Casimir problem then is the usual problem of penetrable, codimension 1 surfaces (see 
[]| for the case with arbitrary codimension). For example, the d = 4 case maps into the IR 3 
Casimir problem with penetrable 2-dimensional surfaces. The limit A — > 00 is the limit of 
perfect conductors (or Dirichlet limit, because = on the conducting surfaces). 

In the rest of the paper we will use this duality to make statements on the nodal line 
statistic from the knowledge of the properties of Casimir energy. 

IV. APPLICATIONS 

Let us start with a well known problem of a Casimir energy calculation: the presence 
of various divergencies, when taking A, A —>■ 00. We will now discuss the interpretation of 
these divergencies for the nodal lines probability. 

A volume divergence oc VA d+1 (divergent when A — > 00) is removed by the factor of 1/Z 
in our definition of P[<p]. This term is however independent of the presence and shape of 7. 
In QFT it would represent a cosmological constant term. For what concerns the divergencies 
arising when A — > 00 (the so-called Dirichlet limit) let us recall that the tube radius r built 
around the reference line 7, goes to when A — > 00. It is then natural that the probability 
S.y(X) of having a nodal line within a distance r from 7 goes to zero when A — > 00. This 
reflects in the fact that £ — > +00 when A — > 00. 
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We also know that the interaction energy (the one that depends on the distance between 
the bodies) remains finite when A, A — > oo. This means that there are some properties of 
the nodal lines, connected with the interaction part of £, which are well-defined also when 
the tube radius r — > and the cutoff goes to infinity. In order to identify them we must 
define a quantity which stays finite in this limit. Led by the intuition about the Casimir 
energy of rigid conductors we recognize that this problem is related to the removal of the 
self-energy for rigid, disconnected bodies. Suppose then our curve 7 is composed of two 
disconnected pieces 7 = 71 U 72 (we always require them to be straight and both parallel in 
the xo direction). Their Casimir energy can be written as 

£7 = £71 4" £72 4" £int (IS) 

where the first two terms are independent on the distance between the nodal lines and the 
last term goes to zero when the distance between the two curves goes to infinity (this can 
be taken as a definition of £- mt ). The terms £ 71 2 are the energies of isolated plates. 

Let us define the quantity V as the ratio of the probability of having ji U 72 and the 
probability of having both 71 and 72 independently of each other: 

Tt '-'71+72 (A) 



e' 1 ^. (19) 



(X)S J2 (X) 

The interpretation of V is the following: if V > 1 (V < 1) then it is easier (more difficult) 
to find a nodal line 72 if another line 71 is present. 

The interaction energy £i nt is always finite (even when m — > and/or A, A — > 00) and 
hence so is V. Moreover we know from quantum field theory that £ int < and that it 
increases when the nodal lines are pulled apart. Hence we can say that the presence of a 
nodal line 71 makes it easier for another nodal line to be born. Hence in this case nodal 
lines induce other nodal lines in their vicinity. 

The choice of the scalar field spectrum allows us to use all the machinery of QFT (includ- 
ing the Hamiltonian formulation) to calculate the Casimir energy £ mt . Depending on the 
value of A, A, m and the distance between the nodal lines a, we can use a weak or a strong 
coupling approximation for £ int . Since A has dimension t d ~ x (£ is a length scale) the relevant 
dimensionless parameter is e = a\ l ^ d ~ l \ Moreover assuming A 3> m, ((V0) 2 ) ~ A 2 we 
have A ~ l/r d+1 A 2 . By choosing a, r, A we have e 1 for weak coupling and e 3> 1 for the 
strong coupling regime. 



We will now make some explicit sample calculation in these two regimes. 
In the weak coupling regime can use a Feynman diagram expansion [8j (one must use the 
Euclidean cutoff A on the k integrals) for the Casimir energy £ 

£ = \£ 1 + \ 2 £ 2 + ... (20) 

where £\ is given by the tadpole diagram, £2 is given by the 2 legs diagram and so on. 
We will now calculate the first two terms of the series (J2*U|) showing that £\ drops between 
numerator and denominator in V and then calculating the first correction to V, i.e. £2- We 
will also show that £2 < 0, which implies V > 1, in a region of order 1/m around any nodal 
line. 

For the tadpole diagram we have 



Ex = J d^xVix) <0 2 (x)) = <0 2 (O)) J d d - l xV{x) (21) 

or in Fourier space 

f d d ~ x k Tr/I . f d d q 1 

where V(k) is the Fourier transform with respect to the d — 1 spatial dimensions 

V{k) = J d A - x xV[x)e lkx (23) 

so that in d = 2 and with V(x) = 5(x) + 8(x — a) we have 

V{k) = 1 + e lka . (24) 

Since J dxV(x) does not depend on a, the tadpole diagram does not contribute to the 
interaction energy E- mt and hence does not contribute to V. We must then go to the next 
diagram, the one with two legs to find the first non zero correction to £- lxsh . The two- legs 
diagram contribution can be written as 

d d - l k Trn ^ r , M f d d q 1 1 

(2^ 



Sl ~ I 77 Hy^ v{k)v{ k) J p^pjn^^' (25) 



It contains an a-dependent interaction term. To calculate this a-dependent term we can 
send A —>■ 00 (for d = 2 we can take this limit safely) and by means of the usual technology 
for handling Feynman diagrams we find 

1 f 1 , f°° dk 



£2 = -— I dx I — e ika — — (26) 

2* Jo ./-oo27r m 2 + x(l-x)k 2 V ; 
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Performing the integrals gives 



£2 = ~ (1 - $(2v^)) , (27) 



where $ is the error function. Then V can be written, to this order in A, as 

<p = e T^(l-$(2 V ^5)) - (28) 

As we said before V decreases when a increases. Moreover for a ^> 1/m we can do an 
asymptotic expansion for $ finding 

V ~ exp (V— ^e" 4 ^ , (29) 



so P ~ 1 effectively for a ^> 1/m. 

The strong coupling limit has to be tackled with different, non-perturbative techniques. 

The d = 2 case can also be solved exactly for any number n of parallel nodal lines by 
using the techniques in |7j. The resulting exact expression for n > 3 is too cumbersome to 



be presented here and we refer the reader to [l] for details. 

Two parallel nodal lines separated by a distance a in the limit A — > 00 are dual to the 
problem of two points in 1 space dimension at a distance a. The Casimir energy for this 
configuration is: 

1 f°° dE ( e -2aVE+W? \ 

£ '-rJ vf^- oTFfT^ j- (30) 

We can use this formula to make some predictions about V. To begin we know that 
£i n t < and that it has a minimum at a = as £mt(« = 0, A ^> A) ~ —A log(2)/27r (here we 
put for simplicity m — 0). It can be proved that this is also equal to £ S i ng i e (2A) — 2£ sing i e (A) 
(where £ S i ng ie is the energy of a single delta function when A — > 00), which appeals to intuition 
since at a = we are just superposing two delta functions to create a delta function with 
double strength. If m > it can be proved that £ oc exp(— 2ma) and hence again V — 1 
when a>l/m. In any case we can say that V decreases when a increases. 

A difficulty must be noticed here, concerning how far one can push the interpretation 
of S 1 as the probability of having a nodal line 7. Reasoning like in |9j, assuming Dirichlet 
boundary conditions on a line 7 intersecting the x\ axis at say X\ = we have then to 
expand our RW in series of sin(kjXi) (with random coefficients). Reasonably the subset of 
RW that has a nodal line on 7 should be expandible in this basis as well. If moreover our 
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spectrum is cut off at A then one expects that for a <C tc/A one should find much fewer nodal 
lines (the first zero of sinz is at z = tt). In fact a similar phenomenon is found in [l| for 
the monochromatic spectrum. The nodal line length density normalized to its asymptotic 
value goes to ~ 0.5 for x\ = 0. However increasing X\ the nodal line length density suddenly 
increases to a value higher than the asymptotic value and then relaxes, oscillating, to 1. In 
analogy our quantity V should then start from a value < 1 at a = 0, increase in a region 
1/A to a value V > 1 and then relax to V — 1. Evidently the first, O (1/A) region is not 
captured by our analysis, while the second one is. This, as we said in the discussion after 
Eq. (|5ji . can possibly be traced back to the failure of the 'mean field' approximation that was 
used to link S" 7 with the true probability of finding a nodal line 4|. We have hence learned 
that we must assume a>l/A for our results to hold. Equation (J30|) for a ^> 1/A, 1/A, and 
m = gives 

£ " = ~Wa < 31) 

yielding 

V = e nT/24a . (32) 

The higher dimension (d > 2) case cannot be solved in general, due to its strong ge- 
ometry dependence. The constant time section of 7 can be any hypersurface representing 
disconnected conductors in space. The Casimir problem is the most generic one and we do 
not posses an efficient way of solving it. We know however how to solve the case of parallel, 
large (actually, infinite) d — 2 hyper-planes (lines in d = 3, planes in d = 4 etc.). The result 
for m = 0, A — > 00 is 

S 

£int OC —r (33) 

ar 1 

where S is the d — 2 dimensional area of the hyper-planes, a their separation and the 
proportionality constant depends on d. 

One of the main problems of Casimir physics is to find effective (analytical or numerical) 
ways of calculating the Casimir energy for arbitrary configurations of perfect conductors. 
Despite recent developments this problem escapes analytical solution for all but the 

simple parallel plates case. We then expect to gain some insights from the other side of the 
duality, namely the nodal lines distributions. 
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V. EXTENSIONS AND FURTHER DEVELOPMENTS 



Extension to different spectra. It would be interesting to know how much of what we 
said in this paper, based on the scalar field spectrum, is valid for other spectra (like the 
monochromatic spectrum). The monochromatic, as well as other kinds of isotropic spectra 
cannot be modelled by an Hamiltonian field theory, even thought their probability functional 
is gaussian. The high degree of non-locality of these spectra implies that the free energy F 
is not extensive. Hence we could not even define a T-independent quantity like the Casimir 
energy S. It is hence of great interest for a field theorist to grasp some of the properties of 
these generalized free QFTs in terms of some, more intuitive perhaps, statistical properties 
of random waves. 

Extension to codimension > 1. Generically nodal lines of real fields have codimension 1 
(lines in the plane, etc.) because they are defined by a single condition, namely <p(x) = 0. 
Codimension 2 or higher nodal lines are non-generic and have extremely low probability of 
occurring. For example the probability of having <p{x) = at an isolated point requires both 
4>(x) = and |V</>(x)| = at the same point. This is extremely unlikely in the sense that it 
las measure 0, and would never show up in a Montecarlo simulation. We know in fact from 
?| that conductors of codimension 2 and higher cannot be defined with A > 0. They must 
be defined as a limit A — > 0~. However the generating functional 5* 7 (A) is not well-defined 
for A < 0. It diverges badly. Actually, since S 7 (X) = (e~ xx ~*}, for A < it is finite and only 
if the probability distribution of A 7 decays at infinity faster than e XXj . It turns out that 
one can take the limits (A — > 0~ and shrinking 7 to codimension > 1) in such a way that 
this divergence and the infinitesimal probability of a codimension > 1 nodal line occurring 
compensate, giving a finite value for S T 

Extension to complex fields and phase singularities. A nodal line of a complex field 
is a more interesting object than that of a real field Complex field nodal lines are 
phase singularities whose strength can be interpreted as a topological charge 5j. Various 
correlation functions of this charges have been calculated by means of the Gaussian field 
technology. It would be interesting to see what the Casimir energy analogy has to say on 
these objects. 

Numerics. One of the main reasons this duality is interesting is that it could lead to a 
more efficient numerical algorithm for computing Casimir energies of conductors of arbitrary 
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shape. However this issue is beyond the scope of this paper and we leave them for future 
work. 

VI. CONCLUSIONS 

We have shown that there is a dual description of random waves in terms of quantum 
field theory. In particular we put forward and started the exploration of the duality between 
the probability of having a nodal line close to a given disconnected reference curve and the 
Casimir energy of a configuration of conductors. 

We used this duality to infer some properties of the distribution of nodal lines and we 
proved that, for the scalar field spectrum, nodal lines induce other nodal lines in their 
proximity. This last statement just follows from the attractive nature of Casimir interactions. 

This duality can be used in the other direction to gain information on the Casimir energy 
of an arbitrary configuration of conductors from the statistical properties of the nodal line. 
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